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We study the evolution of gravitational waves through the preheating era that follows inflation. 
The oscillating inflaton drives parametric resonant growth of scalar field fluctuations, and although 
super-Hubble tensor modes are not strongly amplified, they do carry an imprint of preheating. This 
is clearly seen in the Weyl tensor, which provides a covariant description of gravitational waves. 



O 1. INTRODUCTION 

o : 

^ Gravitational waves in inflationary cosmology are produced by sub-Hubble scale vacuum fluctuations. They are 

stretched to super-Hubble scales by inflationary expansion, and then they generate anisotropies in the cosmic mi- 
crowave background (CMB) (see, e.g., Grishchuk 1975, Mukhanov et al. 1992). Scales which re-enter before matter- 
radiation equality will rapidly redshift (and suffer small damping during recombination), having negligible effect on 
CMB temperature anisotropies, while scales which re-enter later can affect the temperature anisotropies at large an- 
gles. Roughly speaking, we can approximate the impact of gravitational waves upon CMB temperature anisotropies 
by neglecting all local causal dynamics and treating the large-scale fluctuations as frozen beyond the Hubble scale 
until re-entry, where their amplitude is conserved and determined by the inflationary Hubble rate at the time of 
■ Hubble-crossing. This picture is confirmed by the evolution equation (Lifshitz 1946) 
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' for the Fourier modes of the transverse traceless tensor perturbation /y , defined (on a flat background) by 

o 

ds 2 = a 2 (rj) [~dr] 2 + (<% + fe) dx l dx j ] . (2) 

O 

CT^ Equation ([!]) shows that on very large scales, k/aH <C 1 (where H — a' /a 2 ), 

fij « Aij + Bij / — j , (3) 



where A'^ = = B^. If a' > then the B-mode is decaying, and we have /y sw constant, with the constant 
?H ' determined by the Hubble rate at inflationary Hubble-crossing, H(rj c ) = k/a(rj c ). 

This simple picture is modified by small corrections induced at preheating. At the end of slow-roll inflation, the 
inflaton oscillates and transfers its energy to fluctuations, initiating the reheating era, which ends when created 
particles thermalize as a radiative plasma. Reheating is often initiated by a preheating era, marked by coherent 
inflaton oscillations which drive parametric resonant amplification of scalar fluctuations (see, e.g., Traschen and 
Brandenberger 1990, Kofman et al. 1994, 1997). One crucial point about preheating is that, since the inflaton is 
coherent on scales well beyond the Hubble horizon, it is in principle possible for super-Hubble fluctuations to be 
amplified without violating causality (Bassett et al. 1999a, b).n In the case of scalar fluctuations of the metric, this 
can in principle produce nonlinear amplification, depending on initial conditions and coupling strengths (Bassett et 
al. 1999a,b,c, 2000, Ivanov 2000, Jedamzik and Sigl 2000, Liddle et al. 2000). While tensor fluctuations will not be 
strongly amplified by scalar inflaton oscillations, these oscillations nevertheless could leave an imprint on the tensor 
spectrum on large scales. Such a possibility is usually ruled out on causality grounds, but no such causality constraint 
operates during coherent oscillations of the inflaton condensate. 



1 lt is also possible for small-scale gravitational waves to be generated by gravitational bremstrahlung via rescattering of 
scalar fluctuations during preheating (Khlebnikov and Tkachev 1997). This is a particular example of the generation of tensor 
perturbations by scalar perturbations at second order (Matarrese et al. 1998) . 
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These small corrections will be carried into Eq. (|l|) via the scale factor, which inherits an oscillatory addition 
to its average value. However, the nature of the effect is more clearly brought out via an alternative description of 
gravitational waves, based on the idea that a full characterization requires the curvature tensor, not the metric (Pirani 
1957). Transverse traceless modes are given by the divergence- free electric and magnetic parts of the Weyl tensor 

where is the background four- velocity (there are no velocity perturbations for tensor modes). This covariant 
description of gravitational waves was used by Hawking (1966), and is remarkably analogous to electromagnetic 
radiation theory (Dunsby et al. 1997, Maartens and Bassett 1998). 

In this paper, we use the Maxwell- Weyl approach to gravitational waves to investigate the effects of inflaton 
oscillations in some simple preheating models, generalizing previous work in Minkowski spacetimc (Bassett 1997). In 
Sections 2 and 3, we give the basic equations and their qualitative properties. In Section 4 we present the numerical 
calculations, and Section 5 contains concluding remarks and discussion. 



2. BACKGROUND DYNAMICAL EQUATIONS 



The background inflaton is governed by the Klein-Gordon equation 

(p + 3Hcp + V v = , (4) 

where V v = dV/dip. We will consider the simple chaotic inflation potentials V — \m 2 tp 2 and V — jXip 4 . Although 
the resonance in scalar fluctuations is dramatically increased when the inflaton is coupled to other fields, for example 
via the additional potential term ^g 2 f 2 X 2 ^ this does not seem to have an effect on tensor fluctuations (Tilley 2000). 
Thus we will confine ourselves to simplified single-field models of preheating. Preheating in more realistic models ends 
when backreaction effects of the fluctuations destroy the coherence of inflaton oscillations. In our simplified models, 
backreaction effects are not incorporated, but we can use the results from detailed investigations to estimate the time 
that preheating lasts (see, e.g., Kofman et al. 1997). 

The Hubble rate is determined by the Friedmann equation 

H 2 = \n 2 [^ 2 + V(fp)] , (5) 

where n 2 = 8n/M 2 . Equations ([|) and (||) imply H = — ^K 2 Lp 2 . The energy density and effective pressure of the 
inflaton are 

P = k 2 (\£ 2 + V), p = k 2 (i^ 2 - V) . (6) 

During slow-roll inflation, the coupled equations (Q) and (||) have approximate analytic solutions for the simple 
potentials. When the value of if drops low enough, slow- roll ends and the oscillatory regime begins. The approximate 
analytic forms are (see, e.g., Kaiser 1997, Kofman et al. 1997) 

|mV, (?) 
ran) and V = jXp 4 , (8) 

where a,i n = 1 and cn is a Jacobian elliptic function. The initial values of ^„ are ~ 0.3M p in the quadratic case, and 
~ 0.6M p in the quartic case. 

Time-averaging over oscillations shows that a cx t 2 ^ for the quadratic potential, and a oc t 1 / 2 for the quartic 
potential. These are the asymptotic values of the scale factor, but in practice backreaction effects due to couplings 
will end the preheating oscillations. If one uses the average forms a for a, i.e., if one ignores oscillatory behaviour in 
the inflaton, then one regains the standard results for gravitational wave evolution in the dust and radiation eras. In 
particular, models with a smooth transition from inflation to radiation-domination, neglecting reheating dynamics, 
show that there is no super-Hubble amplification of gravitational waves (e.g. Caldwell 1996, Tilley and Maartens 
1998) . Our numerical integrations show that this averaging loses interesting features in the gravitational waves (see 
Section 3). We do not use the approximate forms in Eqs. M) and (H) for our numerical results — instead, we integrate 
the Friedmann and Klein-Gordon equations numerically. 
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3. MAXWELL- WEYL GRAVITATIONAL WAVE EQUATIONS 



Gravitational wave perturbations in the covariant approach are governed by the Maxwell- Weyl equations: 

(divE) M = 0=(div^ , 

+ ZHE^v - curlHfiu = -^(p+p)a fiV , 
+ 3HH^ + curl E^ = , 

where is the shear, a dot denotes it^V^, and div and curl are the covariant spatial divergence and curl for tensors 
(Maartens and Bassett 1998). These equations hold for perfect fluids and minimally coupled scalar fields (in which 
case p + p = K 2 ip 2 )- The shear is a tensor potential for the electric and magnetic Weyl tensors (Maartens and Bassett 
1998): 

Epv = -ct^u - 2Ha flL , , (9) 
H^u = curler^ , (10) 

in close analogy with the Maxwell relations E = —A and H = curl A. Taking the curl and dot of the Maxwell- Weyl 
propagation equations, and using the identity for the tensor curl of the curl, we find wave equations for the three 
tensors (Dunsby et al. 1997, Maartens and Bassett 1998, Challinor 2000). Using equations (||) and (||), these become 

+ 7HE^ + 4k 2 V{v)E^ = AE^ + K 2 <p [5H<p + V'(<p)] , (11) 
+ 7HH^ + 4K 2 V{ip)H^ = AH^ , (12) 

where A is the covariant spatial Laplacian. One can also derive a wave equation for the shear: 

<V + hHa^ + k 2 [2V(tp) - \<p 2 ] = Aa^ . (13) 

We decompose the tensors into modes (Challinor 2000) 

E»u = ^ 2 £ k 2 [£kQ$ + SkQ$] , (14) 
Hp V = a- 2 k2 [HkQ$ + H k Q$] , (15) 
o> = a- 1 J2 k [ S kQ$ + S"Q$\ > ( 16 ) 

where denotes a symbolic sum over harmonic modes, and Q, Q are tensor harmonics of electric and magnetic 
parity, which are time-independent, transverse traceless, and related by 

curlQ« = -Q«, 
p a p 

showing that the different polarization states are coupled. The mode functions £ k determine the tensor contribution 
to the CMB power spectrum. The magnetic Weyl mode functions are algebraically related to the shear mode functions 
by 

H k =S k . (17) 
Specializing the results in Challinor (2000) to the scalar field case, we find the evolution equations 

4 + 3H£ k - | ( ^ - \k 2 ^ 2 \ S k = , (18) 



S k + HS k + -£ k = . (19) 
a 



k 



For comparison, in the time-averaged approximation, the (p 2 term in the coefficient of S k in Eq. ( J18| ) is replaced 
by 7/9, where p — pi„a~ 37 , with 7 = 1 for the averaged quadratic case (i.e. cold matter or 'dust'), and 7 = § for the 
averaged quartic case (i.e. radiation). 
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4. NUMERICAL RESULTS 



Equations ()T^)-(|l9|) for the gravitational wave perturbations, and (|4])-(|5|) for the background are integrated numer- 
ically, and the results are compared with those for the time-averaged approximation. We investigated the evolution of 
COBE modes, which left the Hubble radius at about 50 to 60 e- folds before the end of inflation, so that k/aH ~ 10~ 
at the start of preheating, t = ti n . We also integrated for a typical small scale mode, with k/aH ~ 10 at t — ti n . The 
initial values for ip were given in the previous section. We used the slow-roll relation to set the initial inflaton velocity, 
(pin = — M P V V /V24ttV . For the tensor modes, we used the initial value 10 -5 . The inflaton mass and self-coupling 
were chosen as m — 10 _6 M p and A = 10~ 12 . In order to take account of backreaction effects, which will end the 
coherent inflaton oscillations, we terminated the integrations after t = 100. 

The results are shown in Figs. 1-8. 
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FIG. 1. The electric Weyl mode £k(i~) on super-Hubble scales, k/{aH)i„ = 10 24 . The solid line is for the inflaton potential 
V = \rr?ip 2 , while the dotted line is for the corresponding time-averaged model, which behaves like cold matter (dust). 
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FIG. 2. The magnetic Weyl mode Hk(r), as for Fig. 1. 
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FIG. 3. The electric Weyl mode £k(r) on sub-Hubble scales, k/(aH) in — 10. The solid line is for the inflaton potential 
V — imV 2 , while the dotted line is for the corresponding time-averaged model, which behaves like cold matter (dust). 
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FIG. 4. The magnetic Weyl mode 74 (t), as for Fig. 3. 
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FIG. 5. The electric Weyl mode £k(r) on super-Hubble scales, k/(aH) in = 10~ 24 . The solid line is for the inflaton potential 
V = j\ 4 <p 4 , while the dotted line is for the corresponding time-averaged model, which behaves like radiation. 
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FIG. 6. The magnetic Weyl mode Hk(r), as for Fig. 5. 
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FIG. 7. The electric Weyl mode £k(r) on sub-Hubble scales, k/(aH) in — 5. The solid line is for the inflaton potential 
V = j\ 4 <p 4 , while the dotted line is for the corresponding time-averaged model, which behaves like radiation. 
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FIG. 8. The magnetic Weyl mode Hk(r), as for Fig. 7. 
5. CONCLUSIONS 

Apart from specific features peculiar to the potential, these figures reflect two key properties: (1) Super-Hubble 
gravitational waves are not significantly amplified by coherent inflaton oscillations during preheating; (2) gravitational 
waves on all scales carry some imprint of the coherent oscillatory dynamics of the inflaton during preheating. The 
latter point is clearly brought out by a comparison with the evolution in the absence of oscillations, i.e. for the 
time-averaged scale factor. 

In particular, one can see that the electric Weyl modes on super-Hubble scales, which determine the effect of 
gravitational waves on the CMB (Challinor 2000), inherit oscillations from the inflaton. In the early part of preheating, 
there is also some amplification on average of £k , so that if backreaction takes effect early, then this does produce a 
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small amplification relative to the no-oscillation time-averaged model. However, one cannot expect these preheating 
imprints on £ & to lead to detectable differences in the CMB power spectrum. Firstly, the oscillations will effectively 
be averaged out, and secondly, any amplification is likely to be scale invariant for all measurable anisotropics, since 
the Hubble scale at preheating corresponds to about 1 metre today, so that all cosmologically significant scales at 
preheating effectively have k — 0, and behave like the particular scale chosen in our numerical integrations (Bassctt 
et al. 2000, Jedamzik and Sigl 2000). 

Gravitational waves on sub-Hubble scales oscillate even in the case of time-averaged background, but Figs. 3, 4, 7 
and 8 show that the frequency and amplitude of oscillation are significantly modulated by the inflaton oscillations. 
Thus preheating leaves an imprint on these scales. In principle, this could be detected, but in practice the signal will 
be far too weak, since there is no preheating amplification on these scales. 

The absence of amplification is fundamentally due to the expansion of the universe, since there is strong amplifica- 
tion on a Minkowski background of gravitational waves during preheating (Bassett 1997, Tilley 2000). 
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